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O : Abstract 

(N 

Considering the coherent nonhnear dynamics between two weakly hnked optical waveguide ar- 
n , ' rays, we find the first example of coexistence of Josephson oscillations with a novel self-trapping 

_G \ regime. This macroscopic bistability is explained by proving analytically the simultaneous exis- 

tence of symmetric, antisymmetric and asymmetric stationary solutions of the associated Gross- 
^SJ ■ Pitaevskii equation. The effect is, moreover, illustrated and confirmed by numerical simulations. 

lO ■ This property allows to conceive an optical switch based on the variation of the refractive index of 

^^ \ the linking central waveguide. 
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Introduction. Since its prediction [1| in 1962 and the immediate experimental verifi- 
cation j^, the Josephson effect has led to many implementations in various branches of 
physics. This macroscopic quantum tunneling effect, originally discovered in superconduct- 
ing junctions, is caused by the global phase coherence between electrons in the different 
layers. Similar Josephson oscillatiotis liave been discovered in liquid Helium yl Wi and in 
double layer quantum Hall systems 



The 
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irst realization of a bosonic Josephson junction has been recently experimentally re- 







for a Bose-Einstein condensate embedded in a macroscopic double well potential. 
The difference of the latter from the ordinary Josephson junction behavior is that the oscil- 
lations of population imbalance are suppressed for high imbalance values and a self-trapping 
regime emerges [8|, |9|. The optical realization of a bosonic junction had been theoretically 
proposed much earlier by Jensen [lO| who considered light power oscillations in two coupled 
nonlinear waveguides. 

In order to describe the macroscopic tunneling effect in bosonic junctions, the nonlinear 
dynamics is usually mapped to a simpler system characterized by two degrees of freedom 
(population imbalance and phase difference) while the nonlinear properties of the wave func- 
tion within the single well are neglected. In this approach the symmetric and antisymmetric 



stationary solutions of the Gross-Pitaevskii equation 



wave function 



111 ] are used as a basis to build a global 



13l |. In the weakly nonlinear limit this description is fully valid because 
those symmetric and antisymmetric functions are the only solutions of the equation. For 
higher nonlinearities an asymmetric stationary solution appears, which represents a novel 
self-trapping state ll4J . 

In this Letter, we consider the two weakly linked optical waveguide arrays represented in 
Fig. H When light is mainly injected in one array (e.g. the right one, as shown in Fig^, 
we find a wide light intensity range where it can either remain trapped in this array, or 
swing periodically from right to left and back. The switching from one state to the other 
is determined by a slight variation of the refractive index of the central linking waveguide 
denoted with the index 0. The coexistence of the oscillatory and self-trapping regimes 
corresponds to the simultaneous presence of Josephson oscillation states and asymmetric 
solutions of the Gross-Pitaevskii equation. This is in contrast with all known behaviors of 
bosonic Josephson junctions, where the existence of oscillatory and self-trapping regimes is 
uniquely determined by the parameters of the system. Our theoretical result is expected to 



have a straightforward experimental reahzation in waveguide arrays and in Bose-Einstein 
condensates with double well potential. 
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FIG. 1: Scheme of the proposed setup of two weakly linked waveguide arrays. The refractive index 
of the central waveguide is smaller than the indices of the waveguides in the two arrays. The inset 
displays the elementary cell with the coupling constant Q. 



The model. Waveguide arrays are particularly convenient for a direct observation of 
nonlinear effects, because the longitudinal dimension z plays the role of time. With an 
intensity-dependent refractive index (optical Kerr effect), waveguide arrays become soliton 
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generators, as experimentally demonstrated in [15, 

The model of an array of adjacent waveguides coupled by radiation power exchange is the 
Discrete NonLinear Schrodinger equation (DNLS) 22, |23| and nonlinearity then manifests 
itself by a self-modulation of the input signal (injected radiation). It reads 



i—^ + -nji)j + Q[i)j+i + ipj_i - 2ipj) + li^jl'^ipj = 0, 



(1) 



where waveguides discrete positions are labelled by the index j {—N < j < N), and complex 
light amplitudes ipj are normalized such as to fix the onsite nonlinearity to the unitary value. 
The linear refractive index rij is set to n for all j ^ 0, and to no < n for j = 0. The coupling 
constant between two adjacent waveguides is Q and u and c are the frequency and the 
velocity of the injected light, respectively. Last, we assume vanishing boundary conditions 
ipN+i = "ip-N-i = in order to mimick a strongly evanescent field outside the waveguides. 



The above equation, written for two waveguides (elementary cell in the inset of FigC]), 
straightforwardly reduces to the one considered by Jensen lOll. In that case the resulting 
dynamics is given in terms of Jacobi elliptic functions [9|, |lO| and can be described as 
follows. For a beam of small intensity, light tunnels from one waveguide to the other and 
then back, inducing Josephson oscillations. Indeed, the elementary cell is similar to a single 
bosonic Josephson junction as demonstrated in [8|, |9| and experimentally observed in Bose- 
Einstein condensates J7|. Increasing the injected beam intensity above some critical value, 
light becomes self-trapped and does not tunnel to the other waveguide, which constitutes 
the difference between a bosonic Josephson junction and its superconducting analogue. 

Numerical simulations. We demonstrate now by numerical simulations of model (P) that 
for the device of Fig^ the two regimes, namely Josephson oscillations and self-trapping, 
coexist for a given injected beam intensity and given parameter values. The switch from one 
state to the other is achieved, for instance, by a tiny local variation of the refractive index 
of the central waveguide. 

Let us choose the following values for the parameters in equation (^: 

N = U, Q = 16, Vo = uj{n - no)/2c = 10, (2) 

together with the following input light envelope 

V^,(0) = 0.4sin[(j-iV-l)/5.5], j = l,---iV, (3) 

V,(0) = 0.2 sin[(j + N+ l)/5.5], j = -iV, ■ ■ -0, 

which represents a beam mostly sent into the right waveguide array. Figure |21 displays the 
result of our numerical simulations. While the relative refractive index Vq of the central 
waveguide is kept constant, the power injected initially into the right part of the array 
remains self-trapped, as shown in FigElup to z = 150. A local variation of Vq at 2; = 150, as 
drawn in the inset of Fig|21 makes the self-trapping state bifurcate to a regime of Josephson 
oscillations, which then remains stable. 

Theory. We shall now interpret these results in terms of the continuum limit of 
model (0), which, after an appropriate phase shift, reads as the Gross-Pitaevskii equation 

^^ + ^-V{x)^ + m'ij = 0. (4) 

The wavefunction ip{z, x) depends on the spatial continuous variables z and x = jj \fQ. V[x) 
is a double well potential with a width 1L = {2N + 2)/a/Q and is represented in Fig. 01 The 




Waveguide Number [j] 

FIG. 2: Numerical simulation of the DNLS equation (^ with z = initial conditions ^ and 
parameters 0. By a slight local variation at z = 150 of the refractive index of the central 
waveguide (represented in the inset in terms of the relative barrier height Vq), the regime switches 
from self-trapping to Josephson oscillations. The injected total power is P^ = ^ IV'jP = ^■'^■ 

potential barrier of the double well potential has height Vq and width 21 = 1/^/Q- In the 
region of the barrier we assume, for technical simplification, that the Schrodinger equation 
can be treated in the linear approximation, while numerical simulations are performed with a 
fully nonlinear array. The stationary solution of (JH) are sought as ip{z, x) = $(a;) exp(— z/5z), 
where the real- valued function ^{x) is found in terms of Jacobi elliptic functions as fol- 
lows [241: 



-I < X < -L : $ = 5 cn[yB{x + L) - K{kB), ks] 
-I < X <l : $ = ae^"" + 6e~^^, 
I < X < L : ^ = A cn[7A(x - L) + K(/i;^), A;^], 



(5) 
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FIG. 3: Plot of the double well potential for the continuous model Q): 2L is the well width, Vq 
and 2/ are barrier height and width. The curves are the plots of different types of solutions for the 
total power Pt = 1.7. The inset shows the form of the asymmetric solution for different values of 
the total power. 



with the constants 



lA 



VA^ + P^ IB = VB^ + /3, a = ,/V-f3, 



(6) 



IK denotes the complete elliptic integral of the first kind and the moduli are 

'^A — nf A2 , a\' '^B 



(7) 



2(^2 + /?)' " 2{B^ + (3)' 

The above solutions are given in terms of five parameters {A, B, a, b, (3). As the continuity 
conditions at points x = ±/ provide four relations, the total injected power Pt = J \ip\'^dx 
completely determines the solutions. 

In the weakly nonlinear limit P^ ^ 1, the solutions are either symmetric or antisymmetric. 
The even solution $+(a;) corresponds to A = B in ^ when the solution in the barrier 



region is 2acosh(Ax). The odd solution $_(x) corresponds to A = —B with central solution 
2asinh(Ax). For higher powers, namely above a threshold value, an asymmetric solution 
$a(x) also exists for which A ^ ±B. These analytical solutions are represented in Fig. |21 

Using the symmetric and antisymmetric basic solutions, one can build a variational anzatz 
following 12,ll3| by seeking the solution ip(z,x) under the form 



1p{z,x) = 1pi{z)^i{x) + 1p2{z)^2ix), 



(8) 
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The functions iV'iP aiid |^2p are interpreted as the probabilities to find the light localized 
on the left or right array. By construction, the overlap of $1 with $2 is negligible, namely 



^i^2dx < dx ^l 



dx $2, 



(9) 



where the integrals run on x G [—L, L]. Consequently, the projection of the Gross-Pitaevskii 
equation (JD) on $1 and $2 provides the coupled mode equations [8|, |lO| 



.dip 



dz 



i-EiV'i + I)|V'ir^i = r^2, 



(10) 



D 



j^\dx 
J ^fdx' 



n 



dip2 2 

i^ ^2V'2 + ^1^2! ^2 = ri)i, 

with coupling constant r and nonlinearity parameter D defined by 

"" j^ldx 

The linear levels Ei and E2 are given by 

_ / m^n? + VK] dx 

j^ldx 

and turn out to be equal thanks to ((HI). As a consequence they can be absorbed in a common 
phase in (fTUjl . 

An explicit solution of (fTUI) in terms of Jacobi elliptic functions has been found in [lO| 
and used in Bose- Einstein condensates in 9[. It has a simple form when all the power 
is initially injected into one array, say |'?/'i(0)| = 1, |'?/'2(0)| = 0. Solutions have different 
behavior depending on the value of D being below or above the threshold value 4r 



1,2, 



(11) 



D <Ar : iV'iWP 



D>Ar: \i)i{z)\' 



1 + cn(2r2;, 



Ar' 



, ,D Ar, 
l + dn(-.,-) 



(12) 
(13) 



with |'?/'2p = 1 — I'j/'iP- Solution (jT^ describes an oscillation of light intensity between the 
left and the right array (Josephson regime), since |'?/'i| oscillates around the value 0. The 
period of this oscillation is 

T = 2'K{D/Ar)/r. (14) 



For D > Ar solution (jlHj) oscillates around the value 1/2 + 1/4(1 + a/1 — (4r/D)^) with a 
period 

T' = 8K{Ar/D)/D, (15) 

which shows that the system is in the self-trapping regime. Such simple explicit formulas 
provide two important informations on the response of the system to irradiation of one array, 
namely the period of the Josephson oscillations and the threshold (in terms of parameter 
D) above which we expect a bifurcation to a self-trapping regime. 

However, the picture of possible solutions is not yet complete. We discuss here the 
appearance of a further asymmetric solution when the injected power exceeds a smaller 
threshold value. To be definite, we restrict to the parameter values which follow from Q: 
the width of the double well potential is 2L = 7.5, the barrier width is 2/ = 0.25 and its 
height is Vq = 10. We derive the complete set of solutions © and display the dependence 
of their amplitudes on the total power Pt = J \tp\'^dx in the main plot of Fig0] Below 
the smaller threshold value Pt = 0.76 (showed by the vertical bold dashed line) only the 
symmetric and antisymmetric solutions exist and their amplitudes almost superpose. At 
this threshold value a new solution appears which is asymmetric and the amplitudes in the 
two arrays are represented by the upper [U) and lower (L) branch in Fig. 0] The existence 
of this solution is at the basis of a novel self-trapping regime as we explain below. At a 
bigger threshold value (indicated by the vertical dotted line) solution (jH)) changes from the 
Josephson regime to the usual self-trapping one described above. 

In the region of injected power where the asymmetric solution coexists with the symmetric 
and asymmetric stationary solutions one can induce flipping of one to the other by varying 
the height of the barrier as shown in Fig. |21 It is likely that one can induce such flips by 
other methods (e.g. by the variation of the profile of the injected power). 

Going back to Josephson oscillations related to the symmetric and antisymmetric solu- 
tions, a simple striking numerical check is provided by the inset of Fig. EJ where we plot the 
period of the oscillation against the injected power. Numerical data, which are obtained with 
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FIG. 4: Main plot: dependence of the amplitudes (maximum values on x G [—L,L]) of the sym- 
metric and asymmetric solutions © on the total power (we have not reported the amplitude of 
the antisymmetric solution which almost superposes to the one of the symmetric). The inset dis- 
plays the period of the Josephson oscillations for various values of the injected power: the full line 
corresponds to the analytical expression (|14|) . the crosses are results of numerical simulations. 

the fully discrete (and fully nonlinear) model (Q) for parameters in Q, compare extremely 
well with (HH). 

In conclusion a new coherent state in weakly linked waveguide arrays has been discov- 
ered. This coherent state has the property of being bistable and one can easily switch from 
oscillatory to self-trapping regimes and back. This nontrivial behavior may have interesting 
applications in various weakly linked extended systems, such as Bose-Einstein condensates 
or Josephson junctions arrays, which deserve further studies. 
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